Measurement of inclusive D s , D°, and J/ip rates and determination of the 
production fraction in bb events at the T(5S) resonance 
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The inclusive production of D B , D°, and J/i/j mesons is studied using a 1.86 fb -1 data sample col- 
lected on the T(5S) resonance with the Belle detector at the KEKB asymmetric energy e + e~ collider. 
The number of bb events in this T(5S) data sample is determined. We measure the branching frac- 
tions B(T(5S) -> D S X) /2 = (23.6 ± 1.2 ± 3.6)%, B(T(5S) -» D°X)/2 = (53.8 ± 2.0 ± 3.4)%, and 
<B(T(5S) -> J/ipX)/2 = (1.030 ±0.080 ±0.067)%. From the D s and D° inclusive branching frac- 
tions the ratio f s — (18.0 ± 1.3 ± 3.2)% of to the total bb quark pair production at the 
T(5S) energy is obtained in a model-dependent way. 

PACS numbers: 13.25.Gv, 13.25.Hw, 14.40.Gx, 14.40.Nd 



The possibility of studying B s decays at very high 
luminosity e + e~ colliders running at the energy of the 
T(5S) resonance has been discussed in several theoret- 
ical papers P, @, 0]- Studies of the B s meson proper- 
ties at the T(5S), both alone and in comparison with 
those of its lighter cousins B° and B + , may provide im- 
portant insights into the Cabibbo-Kobayashi-Maskawa 
matrix and hadronic structure, as well as sensitivity 
to new physics phenomena [4j. To date, most studies 
of B s have been performed at hadron colliders, where 
high production rates are tempered by limited triggering 
and detection capabilities. As B factories have amply 
demonstrated for the B° and B + , the kinematic cleanli- 
ness of resonant near-threshold exclusive pair production 
(e + e _ — > T(4S) — > BB), combined with a high trigger- 
ing efficiency and the ability to measure neutral particles, 
can open up a complementary realm of sensitivity to new 
phenomena. The T(5S) could play a similar role for B s 
that the T(4S) has played for B. 

To test the experimental feasibility of B s studies in 
T(5S) events, a sample of 1.86 fb^ 1 was collected with the 



Belle detector over 3 days in June 2005. Earlier an Y[5S) 
data sample of -0.1 fb" 1 was taken at CESR 0, B, 
and, more recently, a dataset of 0.42 fb -1 was collected 
by CLEO [8, 9]. We report here the first results obtained 
by Belle, a determination of the number of bb events in 
the T(5S) data sample, a measurement of the inclusive 
rate of T(5S) events to D s , D°, and J/tp and derivation 
of the fraction of bb events containing B s . 

The Belle detector [l(| has operated since 1999 at 
KEKB llj, an asymmetric-energy double storage ring 
designed to collide 8 GeV electrons and 3.5 GeV positrons 
and produce T(4S) mesons with a boost of f3j — 0.425. 
The recent data sample of 1.86 fb -1 was taken at the 
T(5S) energy of —10869 MeV, under exactly the same 
experimental conditions as in T(4S) and continuum runs, 
except that both beam energies were increased by —2.7%, 
keeping the center-of-mass (CM) boost unchanged. An- 
other data sample of 3.67 fb -1 collected at a CM energy 
60 MeV below the T(4S) just before the T(5S) data tak- 
ing is used in this analysis to evaluate continuum contri- 
butions. 
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Only clean decay modes Df — * (f)ir + ((f) — ► K + K~), 
D° — > K tt^ and — > ^ + ^~ are reconstructed. 
Charge-conjugate modes are implicitly included every- 
where in this Letter. The standard track reconstruction 
and particle identification procedures are used [Io| . The 
invariant mass of — > K + K~ candidates is required to 
be within ± 12MeV/c 2 of the nominal mass. For the 
Df — > decay mode, the helicity angle distribution 
is expected to be proportional to cos 2 6^j; therefore, the 



requirement |cos # hc j| > 0.25 is applied. The helicity an- 
gle is defined as the angle between the directions of 
the K + and momenta in the rest frame. 

In the energy region of the T(5S), hadronic events can 
be classified into three physics categories: uu, dd, ss, cc 
continuum events, bb continuum events, and T(5S) 
events. The bb continuum and the T(5S) events always 
produce final states with a pair of B or B s mesons and, 
therefore, cannot be topologically separated. We de- 
fine the bb continuum and T(5S) events collectively as 
bb events, everywhere in this analysis. All bb events are 
expected to hadronizc in one of the following final states: 
BB, BB*, B*B, B*B* , BBir, BB*_tt, B*Btt, B*B*tt, 
BBirir, B°B°, B°B*, B*B° or B*B*. Here B denotes 
a B° or a B + meson and B denotes a B° or a B~ me- 
son. The excited states decay to their ground states via 
B* -> B 1 and B* s -» B° sl Q. 

An energy scan was performed just before the T(5S) 
data taking to find the peak position of the T(5S) res- 
onance. An integrated luminosity of ~30pb -1 was 
collected at five values of e + e" CM energy between 
10825 MeV and 10905 MeV at intervals of 20 MeV. The 
ratio of the number of hadronic events with i?2 > 0.2 to 
the number of Bhabha events is measured as a function 
of the CM energy (Fig. 1). Here, i?2 is the normalized 
second Fox- Wolfram moment [l3| • This ratio is expected 
to have a Breit-Wigner function shape corresponding to 
the T(5S) resonance, above a flat background. Fixing 
the width value to the PDG value T = 110 MeV/c 2 
(l2l |. the mean mass value is found from the fit to be 
M = (10868±6±14) MeV/c 2 , where the first error is sta- 
tistical and the second error is a systematic uncertainty, 
dominated by the variation of background contributions 
with CM energy. This value is in good agreement with 
the PDG value M T(5S) = (10865 ± 8) MeV/c 2 [l2]. The 
fit value obtained above is treated only as a systematic 
check rather than as a measurement, as the energy range 
scanned is small compared to the T(5S) width, and un- 
certainties due to background contributions are not well 
known. Finally, the energy of 10869 MeV was set for 
subsequent T(5S) runs. 

The uu, dd, ss, cc continuum subtraction method is ap- 
plied to obtain the number of bb events in the T(5S) data 
sample: 
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FIG. 1: The ratio of the number of hadronic events with 
R2 > 0.2 to the number of Bhabha events, as a function of 
the e + e~ CM energy. Only statistical errors are shown. The 
curve is the result of the fit described in the text. 



Here iV|| is the number of bb events in the T(5S) data 
sample, and N^ d and iV c h D a n d t are the numbers of hadronic 
events in the T(5S) and continuum data samples, re- 
spectively. A few percent contribution of t + t~, QED, 
77 and beam-gas events partially cancels in Eq. (1), 
and the corresponding small systematic uncertainty is 
included in the full systematic uncertainty. The effi- 
ciency to select a bb event in the T(5S) data sample, 
e 5S = (99 i and the efficiency ratio for contin- 

uum events in the T(5S) and continuum data samples, 
djsVCSt = 1.007±0.003, are obtained from Monte Carlo 



(1) 



(MC) simulation. The hadronic cross section varies with 
the CM energy as 1/E 2 and a corresponding correction is 
applied. The CM energies for the T(5S) and continuum 
data sets are E 5S = 10869 MeV and E cont = 10520 MeV, 
respectively, with a ~ 5 MeV accuracy of the collider ab- 
solute CM energy calibration. The integrated luminosity 
ratio C$<$,/ Ccont — 0.5061 ± 0.0020 is calculated using the 
standard Belle luminosity measurement procedure with 
Bhabha events. The small statistical uncertainty on this 
ratio can be neglected. 

The value of and uncertainty on the luminosity ratio 
is further checked using high-momentum charged tracks, 
Kg mesons, and D° mesons. To compare T(5S) and 
continuum production, normalized momentum distribu- 
tions are used. The normalized momentum of a particle 
h is defined as x(h) = P(h)/ P max (h), where P(h) is the 
measured momentum of that particle, and -P m ax(M is the 
expected value of its momentum if it were produced in 
the process e + e~ — > hh at the same CM energy. Fitting a 
constant to the distributions of the T(5S) and continuum 
dataset ratios x(h)5s/x(h) cont for 0.5 < x(h) < 0.9 and 
applying small corrections due to the difference between 
the T(5S) and continuum final state particle multiplici- 
ties, which were obtained from MC simulation, the ratios 
0.471 ± 0.005, 0.471 ± 0.005, and 0.477 ± 0.005 are de- 
termined for h = n + ,Kg, and D°, respectively. These 
values agree with the energy-corrected luminosity ratio 
of 0.4740 ± 0.0019 obtained from Bhabha event measure- 
ments, where the factor E^ ont /E^ s was applied to correct 
for the hadronic cross-section energy dependence. 
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From Eq. (1) we obtain the number of bb events in the 
T(5S) data sample, N|| = (5.61±0.03 st at±0.29 syst ) x 10 5 . 
The total systematic uncertainty of ~ 5% includes all sys- 
tematic errors on parameters used in Eq. (1). Finally, the 
bb production cross-section at the T(5S) is measured to 
be (0.302 ± 0.015) nb, in good agreement with the CLEO 
value of (0.310 ± 0.052) nb 0. 

The method for inclusive D s (D s = Df) analysis of 
T(5S) data developed in Ref. [8| is applied here. The 
D s signals {D+ 0tt+,0 -* K+K~) in the T(5S) and 
continuum data samples are shown in Fig. 2(a) for the 
normalized D s momentum region x(D s ) < 0.5, where 
a bb contribution is expected. Here and throughout this 
Letter, the continuum distributions are normalized to the 
T(5S) distributions using the energy-corrected luminos- 
ity ratio. To extract the number of D s mesons, the D s 
mass distribution is fitted by a Gaussian to describe the 
signal and a linear function to describe the background. 
The Gaussian width is fixed to the value obtained from 
MC simulation; the Gaussian mean value and normal- 
ization, and the background parameters are allowed to 
float. The same mass fit procedure, but with the Gaus- 
sian mean value fixed to that obtained from the fit of the 
D s signal in the x{D s ) < 0.5 range, is repeated in each 
bin of x(D s ) in order to obtain the D s yield as a func- 
tion of the normalized momentum, x(D s ). The x(D°) 
and x(J/ip) distributions discussed below are obtained 
using the same fit procedures. 

The normalized momentum x(D s ) distributions are 
shown in Fig. 2(b) for the T(5S) and continuum data 
samples. These two distributions agree well in the region 
x(D s ) > 0.5, where bb events cannot contribute. The ex- 
cess of events in the region x(D s ) < 0.5 corresponds to 
inclusive D s production in bb events. 

The fully corrected x(D s ) distribution for bb events 
is obtained, subtracting the continuum contribution and 
applying a bin-by-bin efficiency correction, obtained from 
the MC simulation. Summing over all bins within the 
interval x(D s ) < 0.5 and dividing by the and 4> decay 
branching fractions and by the number of bb events in 
the T(5S) data sample, the inclusive branching fraction 
B(T(5S) -> D s X)/2 = (23.6 ± 1.2 ± 3.6)% is obtained. 
In the calculations the PDG value B(D+ — » <f>ir + ) = 
(4.4 ±0.6)% [13 is used. The T(5S) inclusive branching 
fraction is multiplied by a factor of 1/2 to compare with 
Bf s \ branching fractions. As explained above, continuum 
bb production cannot be separated from T(5S) events 
and therefore continuum bb production is included in the 
Y(5S)^ D S X branching fraction. The latter is therefore 
defined as the average number of D s mesons produced in 
bb events at the T(5S) energy. 

The dominant contributions to the systematic uncer- 
tainty on the branching fraction measurement are the 
uncertainties due to the B(Df — > 07T+) measurement of 
14%, to the number of bb events of 5%, to the track re- 
construction efficiency and particle identification of 4%, 
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FIG. 2: The D s signal in the region x(D B ) < 0.5 (a) and 
the D a normalized momentum x(D s ) (b). The points with 
error bars are the T(5S) data, while the histograms show the 
normalized continuum (here and in Figs. 3 and 4). 

and to the fit procedure of 2%. 

The obtained inclusive branching fraction agrees well 
with the branching fraction £>(Y(5S) — > D s X)/2 = 
(22 .4 ± 2.1 ± 5.0)% obtained by CLEO [8]. The value 
of £>(Y(5S) — > D s X)/2 is significantly larger than the 
branching fraction B{B -> D S X) = (8.7 ± L2)%, which 
we calculate by combining the PDG average [12j with the 
recent CLEO result [1] adjusted to the value B(Df — > 
(f>ir + ) = (4.4 ±0.6)%. The significant increase of D s pro- 
duction at the T(5S) as compared to that at the T(4S) 
indicates a sizable B s production rate. 

The fraction f a of Bs Bs events among all bb events 
at the T(5S) satisfies the following relation: 

B(T(5S) -» D s X)/2 = f s ■ B(B S -» D S X) + 

(1 - /.) ■ B[B - D S X) , (2) 

where B(B S — > D S X) and B(B — » D S X) are the average 
number of D s mesons produced in B s and B decays, re- 
spectively. Using our measurement of B(T(5S) — > D S X), 
the measured value of B(B -> D S X) = (8.7 ±1.2)% 
and the model-dependent estimate B{B S — > 
D S X) = (92 ±11)% [|, we determine f s = (17.9 ± 1.4 ± 
4.1)%. The systematic uncertainty on f s is obtained by 
propagating the systematic uncertainties on the branch- 
ing fractions included in Eq. (2), taking into account the 
correlation induced by B(Df — * 4>n + ). 

The inclusive production of D° mesons (including both 
D° and D°) at the T(5S) is studied applying a proce- 
dure similar to that used for the D s case. As shown in 
Fig. 3(a), large D° signals (D° — > K~ir + ) are seen in 
the T(5S) and continuum data samples. The number of 
D° mesons as a function of x(D°) is shown in Fig. 3(b) 
for both samples. These two distributions agree well in 
the region x(D°) > 0.5. 

After continuum subtraction and efficiency correction, 
the inclusive branching fraction £>(Y(5S) — * D°X)/2 = 
(53.8±2.0±3.4)% is determined. This branching fraction 
is defined as the average number of D° and D° mesons 
produced per bb event. The dominant sources of system- 
atic uncertainties are similar to the D s analysis, except 
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FIG. 3: The D° signal in the region x(D°) < 0.5 (a) and the 
D° normalized momentum x(D°) (b). 
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FIG. 4: The J/ip signal in the region x(J/ip) < 0.5 (a) and 
the J/ip normalized momentum x{J/ip) (b). 



that the PDG value B(D° -> K-ir+) = (3.80 ± 0.07)% 
has much better accuracy than B(Df — > <fiir + ). 

The value obtained for S(T(5S) -> D°X)/2 is lower 
than the PDG value B(B -> D°A) = (64.0 ± 3.0)% 
as expected, if there is sizable B s production at 
the T(5S). The rate of D° production in B s decays 
can be estimated in a way similar to that described 
in Ref. 0]. Assuming that D° mesons are dominantly 
produced from conventional b — > c processes through 
a fragmentation mechanism with additional uu quark 
pair creation, the branching fraction is expected to be 
B(B S -> D°X) = (8 ± 7)%. Using the inclusive D° 
production branching fraction of the T(5S), B, and B s 
decays and replacing D s by D° in Eq. (2), the ratio 

f a = (18.1 ± 3.6 ± 7.5)% of B^B^ events to all bb 
events at the T(5S) is obtained. The systematic er- 
ror is dominated by the systematic uncertainties from 
B(T(5S) -> D°X)/2 and B(B -> D°X) and is slightly 
affected by the uncertainty on the model-dependent as- 
sumption for the B(B S — » D°X) value. Although the 
D s inclusive analysis provides better accuracy on f s , the 
D° inclusive analysis is an independent method, where 
the uncertainty due to the D° decay branching fraction 
is small and the uncertainty on the number of bb events 
dominates. Moreover, the correlation between f a and the 
number of bb events is positive in the D° analysis and is 
negative in the D s analysis. 

Using the f s value obtained in the D° inclusive anal- 
ysis, the inclusive branching fraction B(B S — > D S X) 
can be extracted from Eq. (2). Using the results of the 
present analysis, and taking into account the correla- 
tion between our measurements of £>(Y(5S) — ► D S X) 
and S(T(5S) -» D°X), we obtain B(B S -> D S X) = 
(91 ± 18 ± 41)%, in agreement with expectations within 
large errors. 

The inclusive production of J/ip mesons is studied in 
the decay mode J/tp — > As shown in Fig. 4(a), 

the T(5S) data sample contains a prominent J/ip sig- 
nal, whereas the J/tp signal in continuum is small. The 
x(J/ip) distributions are shown in Fig. 4(b) for the T(5S) 
and continuum data samples. As expected, J/ip produc- 
tion in the continuum region x(J/ip) > 0.5 is smaller 



than in the low x(J/ip) region where B and B s mesons 
can contribute. 

The difference of these distributions in the region 
x(J/ip) < 0.5 is corrected for efficiency to obtain the 
inclusive J/ip spectrum in bb events. U sing the PDG 
value B{J/ip -> M + M~) = (5.88 ±0.10)% 12J, the inclu- 
sive branching fraction S(T(5S) -> J/ipX)/2 = (1.030± 
0.080 ± 0.067)% is obtained. Systematic uncertainties 
are similar to those in the D s analysis and are domi- 
nated by the uncertainty on the number of bb events. 
The T(5S) branching fraction can be compared with 
B(B -> J/ipX) = (1.094 ± 0.032)% [13], because the 
inclusive J/ip production rates in B and B s decays are 
expected to be approximately equal. Conversely, assum- 
ing the equality of inclusive J/ip branching fractions in 
B and B s decays, we can obtain the number of bb events. 
This cross-check can be important for future large statis- 
tics T(5S) measurements. 

In conlusion, the inclusive production of D s , D°, 
and J/ip mesons has been studied in e + e~ colli- 
sions at the T(5S) energy. The precise measure- 
ment of £>(T(5S) — » D S X) and the first measurement 
of B(T(5S) -> D°X) and B(T(5S) -> J/tpX) branching 
fractions are performed. Clear evidence of sizable B s 
production is observed. The fraction of Bl B s events 
among all bb events produced at the T(5S) energy is de- 
termined from the inclusive D s and D° analyses to be 
f a = (17.9 ±1.4 ±4.1)% and f s = (18.1 ± 3.6 ± 7.5)%, 
respectively. Combining these two f s measurements and 
taking into account the anticorrelated systematic un- 
certainty due to the number of bb events, the value 
f s = (18.0 ± 1.3 ± 3.2)% is obtained. This measurement 
agrees with the CLEO value of f a = (16.0 ± 2.6 ± 5.8)% 
[8| , but has a total relative uncertainty that is a factor of 
two smaller. The B s production rate over all bb events 
at the T(5S) is somewhat larger than the fraction of B s 
mesons produced from the b quark in Z — » bb decays at 
LEP, B{b -> B s ) = (10.2 ±0.9)% 12\. The large f s value 
measured here indicates very good potential for future B s 
studies at high luminosity asymmetric-energy B factories 
running at the T(5S) resonance. 
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